Kerr-Schild type initial data for black holes with angular momenta 



Claudia Moreno 

Center for Gravitational Physics and Geometry, Penn State University, University Park, PA 1680$^ 

Dario Nunez 

Center for Gravitational Physics and Geometry, Penn State University, University Park, PA 16802^ and 
Instituto de Ciencias Nucleares - UNAM, A. P. 70-513, 
Ciudad Universitaria, 04510 Mexico, D.F., Mexican 

Olivier Sarbach 
Department of Physics and Astronomy, Louisiana State University, 
202 Nicholson Hall, Baton Rouge, Louisiana 70803~400lf\ 
(Dated: February 7, 2008) 

Generalizing previous work we propose how to superpose spinning black holes in a Kerr-Schild 
initial slice. This superposition satisfies several physically meaningful limits, including the close and 
the far ones. Further we consider the close limit of two black holes with opposite angular momenta 
and explicitly solve the constraint equations in this case. Evolving the resulting initial data with a 
linear code, we compute the radiated energy as a function of the masses and the angular momenta 
Q • of the black holes. 

o 

00 
(N 



o 



X 



PACS numbers: 02.60.Cb, 04.25.Dm, 04.70.-s, 97.60. Lf 



I. INTRODUCTION 



> 

O 

The collision of two black holes is one of the most promising sources for gravitational radiation which the gravita- 
' tional observatory LIGO will be trying to detect jy . 

The detection of gravitational waves is a remarkable project which present challenges in all of their facets, most 
of them related to the fact that the signal generated from the merging of two black holes will arrive to us as a very 
tiny ripple in spacctimc. Although this is fortunate in the sense that we would not want to have great spacetimc 



O disturbances passing all around us, the fact that the signal is so feeble makes essential to the project to be able to 
^ ' theoretically predict and describe with a great deal of accuracy the gravitational wave generated during the process. 
5-H . The problem of two interacting black holes can be split into three stages, namely: The Newtonian or far limit when 
they are already interacting gravitationally, but so far from each other that the dynamics can be studied within the 
, Newtonian description. Then, as they get closer, the spacetime begins to be modified and no simplifying assumptions 
can be used anymore. It is supposed that the gravitational distortions will increase. This period of the evolution is 
known among the community as the "full 3D" part, meaning that the complete Einstein equations have to be evolved. 
Finally, in the last stage the two black holes are suppose to merge. After merging, the resulting black hole is supposed 
to wiggle and giggle, eventually settling down to a stationary black hole. This last part of the process is known as 
the close limit g], and its importance resides in the fact that it can be described as a single perturbed black hole. In 
this case the evolution can be carried out using linearized perturbation theory. For a distorted Schwarzschild or Kerr 
black hole, the corresponding equations can be cast in a wave-like equation for unconstrained and gauge invariant 
scalar quantities, and the evolution is much more economic, numerically speaking, than the full 3D evolution. 

To solve any system of evolution equations one needs initial data. In general relativity (and gauge theories in 
general), the initial data cannot be given freely but has to satisfy some constraint equations. There is a wide literature 
(see H and references therein) on how to solve the constraint equations of general relativity. Most approaches are 
based on the York-Lichnerowicz decomposition [Q in which case the constraint equations are recast in a set of four 
coupled elliptic equations. These equations simplify considerably if one considers conformally flat initial data, and 
several interesting solutions which describe a slice with two black holes (see Q for a review) have been found in this 
context. However, solutions with conformally flat three-metric suffer from the fact that they cannot exactly represent 
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Kerr black holes since the Kerr metric is not known to possess a conformally flat space part in the line element, due 
to the dragging (for a perturbative proof of nonexistence, see ||.) Therefore, this data is likely to contain "junk" 
radiation which, at least when the two black holes are very far from each other or have merged and settled down 
to a stationary black hole, might not reflect a realistic astrophysical scenario]^]. More recent proposals relax the 
requirement for the initial data to be conformally flat 0, ^) and result in new constructions that are able to naturally 
incorporate spinning black holes. 

In the present work, we construct initial data following an approach introduced by Bishop et al. Q and further 
elaborated in |], |To) which is not directly based on the York-Lichnerowicz decomposition. The advantage of this 
approach is that one can easily obtain closed analytic expressions for the solution to the constraint equations in some 
limiting cases, including the close limit approximation. The approach, which is briefly reviewed in Sec. |j, uses a 
Kerr-Schild type of metric in order to construct the initial data. While in H the superposition of two Schwarzschild 



black holes is discussed, we generalize their proposal to spinning black holes in Sec. III. Our superposition has some 



appealing features (described in more detail in Sec. Ill): i) Far limit: When the two black holes arc infinitely far from 
each other, one obtains in the region near each one of them, initial data describing a single Kerr black hole, ii) Close 
limit: When the black holes "sit on the top of each other", one obtains again a single Kerr black hole. 

In Sec. ^ we consider the close limit regime in which two rotating holes are close to each other with their 
singularities enclosed by a common apparent horizon. In this case, the initial slice can be viewed as a distorted 
t = const. Kerr slice in Kerr-Schild coordinates. We explicitly solve the linearized constraint equations for the case 
in which the angular momenta of the two holes are opposite but equal in magnitude. Using a generalized version |Tl|l 
of the Zerilli and Regge- Wheeler equations, we evolve our initial data in Sec. ^ and compute the radiated energy. In 
particular, we discuss the dependence of the energy on the angular momenta of the holes and compare our results to 



other calculations based on Bowen-York data |12|. We conclude with a summary and remarks in Sec. VI and include 
technical details in appendices |X| and [b]. 

II. KERR-SCHILD INITIAL DATA 

Let us briefly review Bishop et aZ.'s proposal |6) to solve the constraint equations of general relativity: The assump- 
tion is that at the initial slice, the metric and its time derivative are of Kerr-Schild type, 

ds 2 = (i 1 ^ + 2Vl^l u )dx f *dx\ (1) 

where i]^ is the flat metric, i] fJiU dx fl dx L ' — dt 2 — dx 2 — dy 2 — dz 2 , V is a function and a covector which is null, 
normalized such that l t = —1 (as a consequence, S lJ lilj = 1 and U has only two independent components). The three 
metric and extrinsic curvature with respect to a slice t =const. are 

gij = kj - 2Vklj , (2) 

Kij = -i(^/ i ) )t + a[2y/ c (W i / J ), c -(W i ) ii -(^) >i ] ! (3) 



where a = l/y/l — 2V is the lapse and where i — x,y, z refer to Cartesian coordinates. The way Bishop et al. propose 
to solve the constraints is to make an ansatz for the (spatial part) of the null vector , to insert the equations (J|) and 
(||) into the Hamilton and momentum constraint equations and to solve the resulting set of four coupled differential 
equations for the unknowns V, V = dtV and /; = dth- In p], a procedure is given for how to solve the nonlinear 
constraint equations for the superposition of two non-spinning black holes. In the close limit approximation (i.e. in 
the limit where the distance between the two black holes is much smaller than the mass of the holes) the constraint 
equations can even be solved analytically 

Below, we first generalize Bishop et aZ.'s ansatz for k to include the superposition of spinning black holes. Then, 
we solve the constraint equations in the close limit approximation. For simplicity, we will only consider the case 
where the angular momentum of the two black holes are opposite but equal in magnitude. In this case, one obtains 
a perturbed Schwarzschild black hole which is easier to study than a perturbed Kerr black hole. 

III. SUPERPOSING TWO KERR BLACK HOLES 

In Cartesian Kerr-Schild coordinates, the Kerr spacetime can be written in the form (Q) with 

V (4) 

/ _ ( _-\ R x+a y Ry—ax __z_\ (z\ 
M I x > R2+a 2 ' R 2 +a 2 ' R I ' \°) 
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where R is a positive function with units of distance given by 



R 2 = 



i (||f|| 2 - a 2 + ^/(\\x\\ 2 - a 2 ) 2 + 4a 2 z 2 ) . (6) 



Here, M and a denote the ADM mass and the rotational parameter, respectively. Note that R is well-defined outside 
the disc ||x|| 2 < a 2 , z = 0, whose boundary is the ring singularity. 

In order to generalize Bishop et al.'s proposal to spinning black holes, we first notice that the spatial part of the 
null covector given in Eq. (|^) can be rewritten as the sum of two covectors: 



where 



= - 



l w + l {2) 



Rx Ry 



R 2 + a 2 '' R 2 + a 2 ' R 



t 2) = -ay-i.-^o 

\R 2 + a 2 R 2 + a 2 

Here and in the following, we identify li with l l = 6 l Hj and denote by I the vector (l x , l y ,l z ). As in the Schwarzschild 
case 0, the vector can be written as the gradient of a potential function, 

where N is a normalization function, and 

$ = — . 
i? 

Now, the vector K 2 \ which vanishes in the non-spinning case, can be written as the cross product of K 1 ), that is, of 
N V$ (with the same normalization function), and a vector proportional to the angular momentum of the black hole: 



jta) = _jv ± x V*, 
i? 

where a — (0, 0, a). If we allow a to have components in arbitrary directions and replace a ■ z by a ■ x in the expressions 
(Q) and (g), we obtain the Kerr solution for general angular momentum J = Ma. 
Thus, the spatial part of the covector can be written as 



f= A (V$ + V x 6) , (7) 



where the potential $ is given by 



and 



* M 
$= R 



? Ma 
o = 



2R 2 



where now 



R 2 



i (||f|| 2 - ||a|| 2 + v/(||£|| 2 -||a|| 2 ) 2 +4(a-x) 2 ) 



Using the linearity of the V operator, we have the necessary basis to propose the following ansatz for the spatial part 
of the null vector, that is, for the description of a Cauchy slice with characteristics that allow a description of two 
black holes. 
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Let the black holes have masses M\,M2, and angular momenta J\ — M\d\, J2 — M2CI2, respectively and their 
centers be at the points for black hole one, and at x 2 , for black hole two. Having written a single Kerr black hole 
in the form (Q), it is not difficult to generalize Bishop et aVs proposal: We propose 

* = ^ + ^, (8) 

(9) 

with 





Rf = - ^\x-x^-\\h^ + ^[\\x-x^-\\~ ai \\^ +4[a i -(x-x i )j I, (10) 

a, = (l-0(d)) Si + '&id) S f , (11) 

where d = \\x\ — X2W is the "distance" between the black hole's centers, i = 1, 2, and $ (d) is an interaction function 
between the angular momenta of each black hole. Its specific form can be left undetermined for the moment, as for 
the present analysis it is enough that it is a smooth function which behaves in the following way: 

limtf(d) -> 1, 
lim 1? (d) -> 0. 

d — >oo 

There are several well known functions with this type of behavior. The specific form of might be related with 
the spin-spin interactions of particle physics, and will be analyzed in future works. 
The final rotational parameter a/ is given by 

af= M 1+ M 2 ' (12) 

which corresponds to the vector addition of the angular momenta, J/ = J\ + J 2 - We will see shortly that Eq. (jl^) is 
a necessary condition for obtaining the correct close limit. 

Finally, we comment on the fact that our ansatz has potential problems at points where fc = V$ + V x b vanishes 
(such that I cannot be normalized at these points). For the superposition of two Schwarzschild black holes, there 
is exactly one such point. For two Kerr black holes with Mi = M2 and a\ — —02, we found that k is singular at 
where the two holes are centered at x\, x% — —x\. In the situation studied in the next section we 
will see that this singularity does not cause any problems since there, we can excise a region which contains the ring 
singularities of the two holes and the above singularity. 

Besides its mathematical simple form, the above ansatz has the following useful limits: 

1. Far limit: Let us consider that the center of the black hole 1 is at the origin, and that the center of the other 



is far from the origin, that is, ||x2|| goes to infinity. Since then i9 (d) — 0, we obtain from Eq.(ll) that 5j = a,. 
So from equation Eq.([Io|), we see that each R\ depends only on the parameters of the respective black hole. 
Furthermore, in the vicinity of black hole 1, R2 is very large compared to R\ so the second terms in Eqs.(||^) 
can be dropped and one obtains the U for one black hole with mass Mi and angular momentum J\ = M\d\. 

2. Close limit: Suppose that x\ — X2 , i.e. the limiting case when the distance d between the two holes goes to 
zero. In this case, i?(cZ) = 1, and &i = a 2 = at, so R\ = i?2- Taking into account the expression for the final 
angular momentum, Eq.(|l^), the potential and the vector b, Eqs.(H||), we see that one obtains one black 
hole with center at x±, mass M — M\ + M 2 , and angular momentum M a = M\ d\ + M2 3 2 . 

3. One hole limit: If the mass of one black hole, say M 2 , vanishes, it can be seen from Eq.(|l^) that a/ = a% thus, 



from Eq.(ll), 5j — a^. Thus the function R\ only depends on d\ and the potential and the vector b, Eqs. 



reduce to the expressions corresponding to black hole 1. 
4. Schwarzschild limit: Finally, our ansatz reduces to the one given by Bishop et al. in the non-spinning case 

d\ = d 2 = 0. 
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Thus, we consider that we are presenting an educated ansatz for initial data with certain attractive properties 
describing the null vector on a space- like slice with two rotating black holes. In order to complete the initial data, 
we have to solve the constraint equations to yield the function V, and the time derivatives of V and k. We discuss 
this below, where we explicitly solve the constraint equations in the close limit approximation, where the distance d 
between the two holes is infinitesimal. For simplicity, we will also assume that the angular momenta of the holes are 
equal in magnitude but opposite, obtaining a non-spinning black hole when d = (see figure |l|). 



IV. MODEL FOR "KERR+KERR=SCHWARZSCHILD 




FIG. 1: Initial data representing two rotating black holes, such that the merger is a non rotating one. 
We expand $ and b in Eqs. (||J)|), in terms of the dimensionless "distance" between the black holes, 



Mi + M 2 ' 



which we consider to be small compared to unity. In the following, we choose the Cartesian coordinates such that 
the two holes are located at x% — (0,0, M 2 e) and a? 2 = (0,0, —Mie), respectively. This "center of mass condition" 
(Mi xi + M2X2 = 0) on the coordinate system simplifies the calculations below since no terms proportional to e will 
appear in the expansions . Similarly, we shall write Si = Af2.se, and Si = — M\se, where e is a unit vector, and s 
a dimensionless parameter equal to Ji/AfiM 2 . As a consequence of this, 0/ is equal to zero (see Eq. (p2h), thus to 
zeroth order in e, we have a Schwarzschild black hole. 
Using this and the Taylor expansion 



0(e) 



l + ^'(0)e+^"(0)e 2 + O(e 3 ), 



for the "spin-spin interaction" function, we eventually arrive at the following expansions for $ and b 

M M Mi M 2 



b 



r 2r 3 
M Mi M 2 



[3 cos 2 6»- l + q 2 (l-(e-x) 2 )] e 2 + C(e 3 ), 
q cos6»ee 2 + C(e 3 ), 



(13) 
(14) 
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where M = Mi + M 2 is the total mass, and where we have introduced r = ||x||, £ = x/r, cos 8 = z/r. Here, q = si?'(0) 
is proportional to the dimcnsionlcss angular momentum parameter s introduced above. The proportionality factor 
depends on the slope of the interpolation function $(g0 at d = 0. Since $(0) = 1 and $(r) — > as r — > 00 a natural 
choice for this factor is i9'(0) = —1 in which case q = — s. Note that for the superposition of non-extremal black holes, 
we have \s\ < 1. 

To summarize, the expansion depends on the parameters: M-l,M 2 (which are supposed to model the masses of 
the black holes), q = —s (the dimensionless angular momentum parameter) and e (the orientation of the angular 
momentum with respect to the common axis of the holes) . 

In order to proceed, it is convenient to expand the above expressions in terms of spherical harmonics. Since the 
quantities we consider are real, we choose the real harmonics (Y , R lm , I lm 
Parameterizing e = (0, sin x, cos x) without loss of generality and using Eq. ( 
the form (see Appendix [A] for more details) 

-l + 0(e 3 ), 



), m = 1, ...,£, defined in Appendix [A|. 
the spatial vector I eventually assumes 



Ia = 



r 

M 1 M 2 



711 R 11 } 



(15) 



where A — 8, <f> denote the angular variables and e B A denotes the Levi-Civita tensor with respect to the standard 
metric on S 2 . The constant coefficients j£ m and Si m are given by 

1 + y (l-3cos 2 X ) 



72-1 



15 q2 sin ( 2x )' 




(16) 



qcosx, 



qsmx- 



A nice feature of our ansatz is that the only contributions with I < 2 are in the even-parity sector with i = 1, and 
this sector is known to contain only gauge modes (see, for instance, Ref. |ll[] for a proof of this statement). Thus, to 
first order in e 2 , we have no change in the mass nor in the (total) angular momentum. It is also interesting to notice 
that the even-parity contributions look exactly as in the non-spinning case (see || ^|), except that we have to replace 
in each sector M1M2 by j2mMiM 2 (and for q = all the 72m 's vanish except for 720 , in which case we recover the 
ansatz in the non-spinning case.) 

In order to solve the constraint equations, we expand 



V = 

V = 
Ia - 



M 



+ e 2 \v 20 (r) Y 2Q + v 2 -i(r) I 21 + v 22 (r) R 22 ] + 0(e 3 ), 



2 [*ao(r) Y 20 + v 2 - 1 (r) I 21 + v 22 (r) R 22 } + 
2 d A \k (r) Y 20 + k 2 -i(r) I 21 + k 22 (r) R 22 



0(e% 



(17) 



2 -B 



e e 



d B [n 20 {r)Y 20 + n 2 ^(r)I 21 } +0(e 3 ), 



where v 2m (r), v 2m (r) and n 2m (r) are unknown functions. Introducing (|l5| ) and ( |l7] ) into the constraint equations, and 
keeping only terms of the order e 2 , we obtain a set of linear inhomogeneous differential equations for these amplitudes, 
where the inhomogeneity is given by the coefficients in the expansion (|l^). These equations split into two sets - one 
set comprising the even-parity amplitudes v 2mi v 2m , k 2rn , and the other set comprising the odd-parity amplitude h 2m . 
Amplitudes belonging to different m's also decouple because the background is spherically symmetric. 

As observed above, the coefficients in the expansion of Ia in terms of spherical harmonics with even parity are 
proportional to r ) 2m MiM 2 . Since we have already solved the constraint equations in the non-spinning case, for which 
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720 = -\/47r/5 while all other 72m's vanish, the solution to the constraint equations is the same as in the non-spinning 
case H except that we have to rescale in each sector M\M 2 by 72m- The result is 

v 2m (x) = -72m ^ fl-- + -| + -T + ^ ^V^)) + T^^Miy/Mfi), (18) 

where a; = r/M, /1 = MiM 2 /M 2 is the dimensionless reduced mass of the system, and C 2m are free constants. The 
factor 67r in front of the Bessel function I4 ensures that the initial data is asymptotically flat. The amplitudes v 2 m 
and k 2m are obtained from Eqs. (6) and (7) of Ref. g], where, again, one has to replace M\M 2 by r )i m M\M'i. 
In the odd parity sector, the constraints yield 

ri2 m (x) = d 2m , 

x 

where C 2m are other free constants. 

At this point, a question that arises is what meaning can be given to the free constants C 2m and C 2m . As discussed 
in H ^| , the constants C 2m can be used in order to fix the position of the apparent horizon (AH) to its background 
value x — 2: To leading order in e, the position of AH is given by the image of the circle x — 2 under the map 

x x — e 2 D(x) — , 
x 

where the deviation function D(x) is given by 

D(x) = (\v 2m {x = 2)- W ) Y 2m ($,ip). 



The induced metric on the AH is 

gAH = AM 2 (l + e 2 D(x = 2)) (<M 2 + sin 2 1? dtp 2 ) . 
Since the corresponding scalar curvature is 

R ^ = m~2 (2 + 4e 2 i?(.x = 2)), 

we see that the constants C 2m describe the deformation of the AH. There is no physical preferred choice, but for 
defmiteness, we will choose C 2m such that D(x = 2) vanishes, i.e. such that the AH is not deformed. 

A geometrical interpretation of the constants C 2m can be given as follows: In the presence of an axial Killing field 
d v (ip £ (0, 2tt)), we can define the quantity 



J M = 77T [ * (dg w A dxn , (19) 
16tt J s 



where S r denotes a spatial sphere of radius r and * is the Hodge dual. For r — * 00 this gives the Komar angular 
momentum. For an axialsymmetric spacetime with Killing horizon at r — rjj, J{th) gives the angular momentum of 
the horizon. For a Schwarzschild black hole, J(r — 2M ) = 0. Although the Komar formula does not make sense if 
d v is not an axial Killing field, we can linearize Eq. (M) around a Schwarzschild black hole, evaluate the resulting 
expression at the AH, and interpret the resulting quantity as the "angular momentum of the AH" |2C|] . For 1 — 2, the 
result is 

J(r = 2M) = e 2 ^(f> 2m (x = 2) cos$Y 2m dn, (20) 

n m J 

where 4> 2m are the gauge invariant amplitudes given in the next section. However, since cost? is proportional to Y 10 , 
all the spherical integrals vanish by virtue of the orthogonality of the spherical harmonics. Nevertheless, it seems 
reasonable to associate the gauge-invariant quantities <f> 2m {x — 2) with a measure for the differential rotation of the 
AH. Since these quantities depend on C 2m (see below), the constants C 2m can be used in order to control the amount 
of this differential rotation. Again, there is no physically preferred choice for the amount of differential rotation. For 
defmiteness, we will choose C 2m such that <p(x = 2) vanishes. 
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V. NUMERICAL EVOLUTION 



Here, we evolve the initial data obtained in the previous section using the linearized Einstein equations. We compute 
the radiated energy at infinity and compare with other results. Since the background metric is Schwarzschild in Kerr- 
Schild coordinates, we can use the generalized versions of the equations of Zerilli and Regge- Wheeler (RW) derived in 
Ref. 1 11 1 in order to perform the evolution. Expressing these equations in Kerr-Schild coordinates, we have, for the 
case of interest £ — 2, 



x ■ 
x - 



x(x + 2) 
2 

x(x + 2) 



x + 2 
4 



6(3 + 6x + 4x 2 +4x 3 ) _ 
x 2 (x + 2)(3 + 2x) 2 ' 



x ■ 



6(x- 1) 
x 2 (x + 2) 



0-0, 



(21) 
(22) 



where "0 and <f> denote the Zerilli and the RW amplitude, respectively, and where a dot and a prime denote derivatives 
with respect to the dimensionless coordinates t/M and x = r/M, respectively. 

In order to evolve these equations, we have to give initial data to ip, ip and </>, <fi. We can construct those from the 
linearized three metric and extrinsic curvature using the formulas given in Sec. Ill of Ref. p"T[ . The three metric and 
extrinsic curvature can be obtained from the expressions (|^, |^) . The final result is 



ip2m = M 

1p2m = ~M 

where V2 m is given in Eq. (U9), and 



x 3 v 2m - 672mA* 



3x(2x + 3) 

2x 3 (2ie - l)v 2m + x 4 (x - 2)v' 2 



2in 



+ &l2m^{x - 1) 



6x 2 (2x + 3) 



-2M/J,S 2 m 

2MnS 2m 



G 



2m 



4X 



Here, we choose the constants C 2m and C 2m to be 

C 2m = 2.7, 



C 2 



(23) 



which insures that the AH is not deformed and shows no differential rotation, as discussed in the previous section. 

An expression for the radiated energy in terms of gauge- invariant quantities is given in p"l) . In our case (£ = 2) 
this energy takes the form 



dE 
du 



2ttM 2 



E(4 



2 

2m 



^2m 



where ip and <f> are evaluated in the radiative zone. (Note that the radiated angular momentum is zero in our case 
since both the initial and the final stage have zero total angular momentum.) Given the linearity of the equations 
and the form of the initial data, the dependence of the time derivatives of the Zerilli and RW functions are 

^2m{t,r) = e 2 iJ,Mj2m i>o(t,r), 
4>2m (t,r) = e 2 [iM 6 2m <fio(t,r), 

where the functions ipo and 4>o are dimensionless. Accordingly, the total radiated energy is 

~fJ> [Ceven E even + C dd E ot id} , 



(24) 



where 



= J>o 2 !f 




odd 







len ~ 2^r (^20 + 72-1 + 722J ' 
Codd = & (5| + ■ 
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Using the expressions (|l6|), we obtain, after some simplifications, 



1 - y ) + I 2 sin 2 x 



15 



[4 - sin 2 X ] 



Since E even and E a dd do not depend on the angle Xj w c can easily derive the expressions for the extrema of the 
radiated energy, Eq.(p4]) with respect to the angle x- 



dE 
dx 
d 2 E 
dx 2 



6f_ 

5 



MeV E even ^(1--T| sin(2 X ), 



Me 4 fi 2 E even fl - 1 T ] < < ,s ( 1 \ ) . 



where T 



J'° dd . Thus, we see that the extrema are at x = 0, ^, and provided that E even < ^E dd, the maximum 
of the energy radiated will be when the angular momentum of the black holes are parallel to the approaching axis, 
and the minimum will be when the angular momentum is perpendicular to it. Fixing Ci m and Ci m according to Eq. 
( p3| ) , a numerical evolution using a first order reformulation of the equations ( pT|j2^ ) described in Appendix [b| yields 
(to 1.5% accuracy) 



E R „ Rn = 1.74 x 10" 
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E odd = 1.78 x 10 



-3 



so the minimum of the radiated energy is at x — 77 / 2- The fact that 4^0^/9 is much larger than E even suggests that 
this will also be the case for different values of Ci m and Ci m , as long as these do not differ much from the values 
chosen here. 




FIG. 2: The radiated energy as a function of the angular momenta for two cases: The upper curve corresponds to the case 
when the angular momentum of the black holes are parallel to the approaching axis (% = 0), and the lower one when it is 
perpendicular to it (x = 7r/2). Here, we have set e = 1 and fi = 1/4. 

Notice that once E even and E dd have been computed numerically, the radiated energy can be obtained for any 
values of the angular momentum number q and the angular momentum orientation angle x by using Eq. ( |24| ) . This 
function is illustrated in the figures |^ and [| The dependence on the angle x is much more noticeable when the 
angular momentum per square mass of the black holes is large. The maximum radiated energy for a fixed total mass 
and separation is given by 



E < E(q = 1, x = 0) = Me 4 x 5.4 x 10" 
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FIG. 3: Dependence of the energy radiated on the angle between the angular momenta of the black holes and the approaching 
axis, for values of q close to the extremal case. Here, e = 1 and fi = 1/4. 

Our results on the energy radiated agree qualitatively with those presented recently by Khanna |j[2| . There, the 
physical situation sketched in figure [j] is modeled using Bowen-York initial data in the close-slow approximation. It is 
interesting to notice that the quantities we have identified with a measure for the differential rotation of the AH also 
vanish in the initial data used in ]l^ since there, perturbations with odd parity are trivial. In particular, Khanna 
also obtains that there is more radiated energy for x — than for \ — ir/2. However in fj^ , the initial data depends 
linearly on q, so the radiated energy depends only quadratically on q whereas in our initial data, the dependence on 
q is more complicated. In particular, we have a contribution in both the even and the odd parity sectors. 

VI. SUMMARY AND CONCLUDING REMARKS 

We have made an ansatz to give initial data representing the superposition of two Kerr black holes. The ansatz 
depends on the parameters Mi (masses), Ji (angular momentum), and an (unphysical) distance d between the two 
black holes. The ansatz is such that when d goes to infinity, one obtains two Kerr black holes in Kerr-Schild coordinates 
with masses Mi and angular momentum Ji. Furthermore, when d = 0, one obtains one Kerr black hole in Kerr-Schild 
coordinates with parameters M = Mi + M 2 , and J = J\ + J 2 - 

The ansatz superposes two black holes with vanishing linear momentum, and a natural question is whether or not 
we can generalize our proposal to boosted black holes. In their paper Bishop et al. indicate how to superpose 
two boosted Schwarzschild black holes: There, one can replace the potential $(x) = \/\x\ by a more complicated 
expression $#(;?) that depends on the speed, v, of the boost. A proposal for two boosted non-spinning black holes 
which is along the lines of this article would be to take 

= Mi$^ (x - xi) + M 2 $~ 2 (x - x 2 ), 

where Vj are defined in a similar way as the dj above. Whether or not this can easily be generalized to spinning black 
holes remains an open question. 

For the case of two un-boosted Kerr black holes which have opposite angular momentum but which are equal 
in magnitude, we explicitly solved the constraint equations in the limit in which the black holes are close to each 
other. When solving the constraint equations, integration constants appear. These are associated with a deformation 
and a differential rotation of the apparent horizon. Although we have a geometrical picture for them, we do not 
know which values these constants will take in an astrophysical realistic collapse. It should be stressed that this is a 
generic problem of close limit calculations since realistic data can only be provided once the full 3D part of a black 
hole collision is solved. Fixing the integration constant to some specific value, we evolved the resulting distorted 
Schwarzschild black hole using the generalized Zcrilli and RW equations in arbitrary coordinates and computed the 
radiated energy. Our result might change if the constants are chosen differently, but nevertheless, when compared 
to results with a different choice of initial slice, we found a qualitatively similar behavior of the radiated energy as a 
function of the angular momentum of the black holes. In particular, we confirmed that there is more energy radiated 
when the angular momentum of the black holes are parallel to the approaching axis. 

A similar analysis to the one presented here, but for the far limit case, is already underway and will be helpful in 
giving us a clearer physical understanding of the parameters appearing in our initial data, in particular in view of 
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post-Newtonian calculations. Also, the physical meaning of our initial data in the intermediate region remains to be 
understood. 

Finally, we plan to repeat our results using a version of the Teukolsky equation in horizon penetrating coordinates 
|P^t . 0] in order to test the latter and in order to proceed with more confidence to the case where the final black hole 
is a rotating one. 
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APPENDIX A: HARMONIC DECOMPOSITION 



When expanding real quantities defined on the two-sphere S 2 , it is convenient to work with the real spherical 
harmonics. In terms of the standard harmonics Y im , these are defined by 

~ V2 [ h 

As the standard harmonics, the real harmonics defined here provide an orthonormal basis of square integrable functions 
on S 2 . 

For I = 1, we have 



R = \ — sin 9 cos < 

4-7T 



I = \ I — sin 9 sine 

47T 



Note that with respect to this, the radial unit vector x can be written as x = y/br/3(R u ,I n ,Y 10 ). For I = 2, the 
real harmonics are 



y20 

R 21 
I 21 
R 22 

j22 



5 1 

47 2 
T5 
4tt 
15 



(3 cos 2 9 - 1) 



sin u cos tt cos ( 



sin tt cos tt sin < 



15 

— — sin 2 6'cos(20), 

1D7T 



— sin 2 6»sin(2^)). 

167T 



Choosing e = (0, sin cos \) i and expanding the expressions for $ and b given in ( 13 lj) in terms of the harmonics 
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above, we obtain 



$ = — 



M MM 1 M 2 



r r 

,2 



l + -L(l-3cos 2 X ) 



20 



1 



[2 sin X cos x/ 21 - sin 2 xi? 22 ] j e 2 + 0(e 3 ), 



b ■ dx = — 



M Mi M 2 4ir 

3 q \ T 

r 6 V o 



cosx 



V3 V 15 1 



2n , smx /2 i 



^/5 



dr 



smx 



IV5 



dl 21 + idR 11 



cos X dY A > > e 2 +0{e 6 ). 



where * denotes the Hodge dual on S 2 . Using the above expressions and rewriting Eq. (Q) as 

T- dx = N (d$ + *d(b- d&)\ , 

where * denotes the Hodge dual on the three dimensional Euclidean space and N is chosen such that I is normalized, 
it is not difficult to find the expressions given in (|l5|). 



APPENDIX B: NUMERICAL EVOLUTION 



Let us consider the following general one dimensional second order differential equation: 



dt u + 2 /i d t d x u + f- 2 d 2 u + f 3 d t u + f 4 d x u + f s u + f 6 = 0, 



(Bl) 



with u and fi functions of (t, x). We want to reformulate this equation as a system which is first order in time and 
space. To this purpose, we introduce the following new variable: 



7r = (d t + g+ d x ) u, 



with g± to be determined below. We have 



{d t + g- d x ) 7T = {d 2 + (.<?+ + g_) d t d x + g+ g_ d 2 + [(g+). t + g- (g+), x ] d x } 



Comparing this last expression with Eq. (Bl), we see that if 

9± 



the original system (Bl) assumes the form 



9 t v = Aa x v + Bv + S, 



with 



where 



-9+ 

A) -g- 



B 





-k 



A 

Applying the above to the equations (|2l] 

.9+ = 



= (g+),t + 9- (g+),x + hg+ 

22), we find, in both cases, 
x-2 



9- 



A = 0. 



(B2) 



Since A is diagonal with real eigenvalues, the evolution system (B2) is strongly hyperbolic. One of the advantages 
of having a strongly hyperbolic system is that there is a well-defined way how to give boundary conditions (see, for 
example, ]17[]): u and ir are characteristic modes with characteristic speeds — g + and — respectively. Inside the 
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horizon (x < 2), both speeds are positive which means that all modes leave the domain if we put the inner boundary 
inside the horizon. Thus, we can excise the singularity and extrapolate u and tt at the inner boundary, which we 
put at x < 2. As long as x > 2, u is a mode with positive speed and it a mode with negative speed. At the outer 
boundary, we impose the outgoing wave condition n = 0. 

We have written a code that numerically evolves (B2). This code is a second order accurate finite differencing 
one. Spatial derivatives are discretized using centered differencing and the time update is performed with an iterated 
Crank-Nicholson algorithm using two iterations. We have tested our code by comparing the resulting waveforms with 
the ones produced by the code used in Ref. ||. We have calculated the quantities E even and E Q dd defined in Sec. |v| 
by integrating the Zerilli and Regge- Wheeler function, respectively, over time at a fixed location r b s far away from 
the event horizon of the Schwarzschild background black hole. We have checked that the numerical results for E even 
and E a( id at fixed r D b s converge as numerical resolution is increased and that these results do not change by more than 
1.5% when r b s is increased (with typical values of r b s lying in the range 101 M < r b s < 151 M.) 



See the site |ittp: / /ligo. caltech.edu , for the most recent advances on the LIGO. 



R. Price and J. Pullin, Phys. Rev. Lett. 72, 3297 (1994); P. Anninos, R. Price, J. Pullin, E. Seidel, and W.-M. Suen, Phys. 
Rev. D 52, 4462 (1995). 
G. Cook, Living Rev. Rel. 3, 5 (2000). 
J.W. York Jr., J. Math. Phys., 14, 456 (1973). 
A. Garat and R. H. Price, Phys. Rev. D 61, 124011 (2000). 

N. T. Bishop, R. Isaacson, M. Maharaj, and J. Winicour, Phys. Rev. D 57, 6113 (1998). 
R. A. Matzner, M. F. Huq and D. Shoemaker, Phys. Rev. D 59, 024015 (1999). 
S. Dain, Phys. Rev. Lett. 87, 121102 (2001). 

O. Sarbach, M. Tiglio, and J. Pullin, Phys. Rev. D 65, 064026 (2002). 

G. Khanna, Phys.Rev. D 65, 124018 (2002). 

O. Sarbach and M. Tiglio, Phys. Rev. D 64, 084016 (2001). 

G. Khanna, Phys. Rev. D 63, 124007 (2001). 

A. Ashtekar, C. Beetle, O. Dreyer, S. Fairhurst, B. Krishnan, J. Lewandowski, and J. Wisniewski, Phys. Rev. Lett. 85, 
3564 (2000); A. Ashtekar, C. Beetle, and J. Lewandowski, Class. Quantum Grav. 19, 1195 (2002); Phys. Rev. D 64, 044016 
(2001). 



A. Ashtekar and B. Krishnana, Dynamical horizons: energy, angular momentum, fluxes and balance laws, gr-qc/020708C . 
C. Moreno and D. Nunez, Int. J. of Mod. Phys. D 11, 1331 (2002). 
C. Moreno, D. Nunez, and O. Sarbach, in preparation. 

H. Kreiss and J. Lorentz, Initial boundary value problems and the Navier Stokes equation, Academic Press, San Diego, 
(1983). 

Here, we assume that the final black hole is a Kerr black hole. 

The radiated energy, being a coordinate invariant quantity, does not depend on this coordinate choice. 
A more precise way to calculate the angular momentum of the AH would be to use the corresponding formulas for an 
isolated horizon or the newly introduced notion of dynamical horizon E3] . After linearizing the expressions in jh| , we 
have obtain the same result as in Eq. (|2C|). 



